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Abstract 



In this article, we assume that a cold charged perfect fluid is constructing a spherical relativistic star. 
Our purpose is the investigation of the dynamical properties of its exterior geometry, through simulating 
O . the geodesic motion of a charged test particle, while moving on the star. 
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> ; 1 Introduction 

Einstein field equation leads to the solutions satisfying the definition of a perfect fluid, a homogenous 
] isotropic mass and energy distribution with no viscosity, such as Friedmann-Lemaitre models (for a good 
T— I ■ review see [Ij). Moreover, It also provides the chance to investigate the interior distributions of the perfect 
. fluid through energy-momentum tensor; not to forget that the energy-momentum tensor can actually be 
I ] extracted from charged bodies bearing interior pressure. 
^ ■ A relativistic star is supposed to consist of such a fluid either charged or with axial rotation j2]. In 

• ^ ■ order to find out the relativistic, isotropic mass and charge the interior solution has been suggested by B. 
Guilfoyl 13J. 

■ In this article, we are about to form a static charged relativistic star using these so-called solutions; 

besides, studying the dynamical properties of the exterior geometry via a rotating test-particle is included. 

As it goes without saying, in general relativity geodesies motions are of great importance. Since, it 
makes precious anticipations about the cases such as precessions of periastron in planetary orbits and 
theoretical explanations for light deflection in gravitational fields. 

Furthermore, there are lots of researches which have been dedicated to the study of geodesic motions 
in different kinds of geometrical background, which could be derived from general relativity, for example, 
Schwarzchild, Schwarzchild-de Sitter, Reissner-Nordstrom and Kerr geometries (see @]-[9] and an important 
text book [To]). The exact equations of motion have been also derived [TT], and lots of further theoretical 
predictions provided (see [l2]-p^). 

In this case we are considering Reissner-Nordstrom (RN) geometry as the exterior geometry of a static 
charged relativistic star. The procedure that we are going through in this paper concerns the process in 
which charged particles get trapped in gravitational and electrical fields of relativistic stars. 
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The chosen approach is a Lagrangian formahsm and numerical simulations also have been involved. In 
order to do so, three classes of solutions presented by Guilfoyl have been utilized to construct a relativistic 
star and to compare the pertinent geodesic behaviors of the test charged particle. 

The article is organized as follows: Section two contains dynamical preliminaries for a charged particle 
in RN geometry, Section three considers a brief review on Guilfoyl's interior solutions. In Section four, we 
introduce the Lagrangian formalism. Finally, Sections five and six, consist of the simulations of charged 
particle's motions on various types of relativistic stars. 



2 Charged particle in RN geometry 

First, let us consider a spherical source which has been schematically illustrated in Figured) It indicates 
the charged star, considered to be the gravitational source, with Rq being the radius, and the total mass 
and charge being mo and go respectively. A test-particle which is moving on this star, is assumed to travel 
on the RN geodesies. The RN geometry is a static spherically symmetric vacuum solution for the Einstein 




Figure 1: A schematic illustration of a charged spherical star and rotating test charged particle, 
equations in for dimensions {d = 4) and is defined by the metric: 

ds^ = -c^A{r)dt'^ + B{r)dr'^ + r'^dO'^ + sin^ 6'#^ (1) 
in which the variable coefficients are: 

- 1 - 1 2Gmo ^ qlG 
B{r) rc^ 47reor-^c^ ' 

where G stands for the gravitational constant and c for the speed of light. In order to simplify the 
calculations, we use geometrical units, fixing: 

G = c=-^ = 1. 
4vreo 

According to these units, one can relate length dimensions to quantities like mass and electric charge, 
i.e. (see [H] appendix 4): 

[mass] = [mo] = 
[electric charge] = [go] = ^, 
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where £ notates the length dimension. Since c and G are supposed to be 1, the remaining ratios 

2 

and , become dimensionless. These assumptions are commonly used in general relativity. Taking these 
into account, one can rewrite ([1]) as: 

= -(l - ^ + + (l - ^ + ^ly'dr' + r'dn\ (2) 

We will use this metric in further considerations. 

If we write the Killing equations for Killing vector , as pO] : 

(V^X), + (V,X)^ = 0, (3) 

one can obtain 10 Killing vectors for the spacetime defined by metric ([2]). Therefore, since for the 
RN spacetime we have (1 = 4, the number of Killing vectors would be '^^'^^^^ . Due to this fact, the RN 
spacetime has Killing symmetry and we can define energy and obtain time- like geodesies. 

By considering the characteristics of the test-particle (m stands for the mass and q for the charge), 
one can obtain the potential and energy conditions, using the Hamilton- Jacobi equation of wave crests |21] : 



g^''{Pf, + qA^){P,+qA,)+m^ = 0. (4) 
Pfj. is the momentum 4-vectoi0 

= g^aP" =9f.a^- (5) 

Here A is the geodesies affine parameter of the trajectory. The metric components g^^, can be derived 

from the exterior geometry of the source, the RN metric ([2]). Also the vector potential for the static 
charged source will be: 

^^ = ((^(r), 0,0,0), (6) 

where ip{r) = ^ is the scalar electrical potential, outside the star, producing a scalar filed within the 
spacetime. One can define the two conserved quantities as: 

E = -Po, (7) 

the energy, and 

L = P^ (L > 0). (8) 

the angular momentum. Here we consider 9 = ^ , for which the particle's motion is confined to the 
equatorial rotations. Therefore 

Substituting the metric (l2|) in (jH yields: 

(^-220 ) 2 2mo qL ,,dr^, , 

^ +1 - + ^ r bv + ^ + = 0, 

i ^ + 

or 

i^f = (i^ - ^)2 - (1 - ^ + 4)K + (9) 



^Here we use the notation in the first of Ref. [21] . in which in §25.3 the 4-momentum has been defined like Eq. ([5]). 
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Equation ([9]) can be rewritten as: 



[E 



-mo 



2mo 



L2, 



We introduce the potential, felt by the test-particle as: 

990 



V{r) 



2mo g5 



+ 



(10) 



to be assured that positive potential is available. To discuss the motion of the test-particle also, we need 
the rate of variation of rotation angle (p with respect to r . Previously we defined: 



Substituting in Eq. Q yields: 



. dr , 



L2 



= L 



mo. 2 

r 



d\ 



L 



(1 



(11) 



solutions to equations like ()lip . have been obtained using the Weierstrass function and elliptic integrals in 
however in this context, we use a Lagrangian formalism instead of the geodesic equations. 

For different values of L and E , the test-particle experiences different types of motion. The motion of 
the particle, depends on its energy, for sets of solutions for E — V{r) = . The RN spacetime allows three 
types of orbits; for one zero, we will have scape orbits, for two zeros, periodic bound orbits are available, 
and for three zeros, we will have scape/capture orbits and periodic bound orbits. 

It is also worthy to investigate the potential that the particle is going through. The potential may have 
some extremum points, as below: 

dV{r) 

We can rewrite the potential in Eq. (jlOp as: 

V{r) = C{r) + 
and g{r) = 1 



0. 



(12) 



g{r){m? + ^). 



We have C(r 
equation: 



220 
r 



2mo _|_ ql 




(13) 

Using (fT3]) in (fT2]) forms the following differential 



m 



0. 



Solving the above equation for L yields: 



L 



2 r^gC'^ - r^g'^m^ + 2g' m^gr ± 2 V^'^ - 2 r^g^C'^g' + A g^r^m'^C'^ 

r'^g''^ — Ag' rg + Ag"^ 



(14) 



This solution, corresponds to the stable orbits of the test-particle around the star, where the potential 
has its extremum points. The expression for L in Eq. (jl4p indicates that, for stable orbits regardless of 
the energy of the particle, the total angular momentum, for definite values of mass, charge and radius of 
the source, can opt only definite values. As we will see, the energy of the particle is peculiar, when we 
work with a particular potential. 

We are going to discuss the types of motions in Section 5, considering different interior mass and charge 
distribution for the star. We shall introduce these interior solutions in the next section. 
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3 Weyl-type interior solutions for a charged perfect fluid 

In this article, a relativistic star, is defined to be a charged sphericahy symmetric cold perfect fluid. Such 
a fluid is a solution of the Einstein field equations 

R„u-^g^^uR = 8^^{T„^ + E^,u). (15) 

and the Maxwell equations 

VuF^"" = 47rJ'^, (16) 

for a charged isotropic mass distribution. In Eq. (jlSp . E^j^^, is the electromagnetic stress-energy tensor, 
defined as [2T] : 

AttE^, = FPF.p - -^g^.Fp^FP'' , (17) 

in which 

is the antisymmetric electromagnetic field strength tensor. The energy-momentum tensor for an 
isotropic mass distribution with pressure is: 



T 



/p 0\ 

p 

p 

\0 pj 



(18) 



where p is the matter density and p is the fluid pressure. The following spherically symmetric metric is 
assumed to explain the interior geometry of the perfect fluid: 

ds^ = -u^dt^ + ^{rfdr^ + r^dn^, (19) 

in which the metric potential 

a;((^)2 = a + 6^9 + (^^ (20) 
by Weyl's assumption [23], only depends on the electric potential (/?(r) , and a and b in (pO|) are constants. 



The non-zero components of Maxwell equations (|16p for metric (I19p . determine the total charge inside 
a volume with radius r . We have 



Q{r) 



An [ a{r)r'^C{r)dr, (21) 
Jo 



in which cr(r) is the volume charge density. Since outside the star, the geometry obeys the RN metric 
([2]), Equation (pT]) can be rewritten as 



Q(r) = r'^ip'{r)[A{r)B{r)]-2. 
According to the Weyl's assumption in ([20]) . this gives: 

^yb^ + 4A{r) - 4a 

Therefore 

^2 



Q{r) = T^'(r){A(r)S(r)(62 + 4A{r) - 4a)}"i (22) 
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Now to realize what the interior solutions of a relativistic star really are, let us consider another star, 
with radius tq {vq = Rq + 6r), charge Qq and mass M which are considered to be known values. This 
star (Figure ED is not the one that will be used as the massive source for the test particle to move on. It is 
only used for discussing the interior solutions and it is in fact, a boundary condition to derive the desired 
constants in the interior solutions. However, the massive source in Figure [T] has the same distributions but 
in lower dimensions. The radius rg will be considered to be the same for all cases, therefore the stars under 



Qo,M 




Figure 2: A star with higher dimensions, having the same distributions as Figure [U but with known values for 
charge and mass. 

discussion, will have the same geometrical appearance and obviously the amounts Qo and M will differ 
for different interior solutions. Now for the star in Figure El the continuity condition at r = tq implies 
that: 

B[ro) ro 
in which Qq = Q(r'o) . Therefore, from Eq. (|22p we obtain: 

M 1 



±^^62 + 4(l-aQ2), (23) 
which implies that for a = b = we have: 

M = \Qo\. 

This is the condition concerning a Majumdar-Papapetrou star for a Weyl-type metric potential. The 
pressure for these stars is zero; in this case the relation between the metric potential in (j20p and 
the electrical potential (p{r) is always a perfect square \25\ I26j . Also we will discuss the motion of the 
test-particle around these stars in Section 5. 

In Ref. [3], Guilfoyl has presented some interesting Weyl-type interior solutions for a spherically 
symmetric charged perfect fluid. He based his solutions on two presuppositions: 

• The Schwarzchild condition which implies: 

3 

SttToo — 87r£^oo = — n = constant, 

or 

87rp(r) + ^ / = TPJ = constant. (24) 
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• The expression for the radial part of metric (jl9p : 

= (1 - (25) 

In these two presuppositions, the quantity , is a constant that must be determined using the bound- 
ary conditions at r = vq (in Figure [2|) . We now bring up two classes of Guilfoyl's solutions which we will 
use later to form the test-particle's potential |3j: 

3.1 The case of Class I 

In this case, the metric potential is considered to have the spherical form: 

au^ = b + {c + ipf, (26) 

where c is a constant. We call this, the Weyl-Guilfoyl metric potential as it is in [24J. Both using the RN 
limit in (|19p and also considering (j26|) we have: 



As we can see, while a = 1 and 5 = we have M = \Qq\ which is the Majumdar condition. For the 
expression (j26p . two sets of solutions have been introduced: 



Class la : a = | , b = 



The function ^(r) is defined as: 



Q{r) 



. kr'^ 



(28) 



, 1 



^l;(r)=l-kR'{l-j^) 



(29) 



Also the constants , k and / can be determined, using the boundary conditions at r = tq (Figure 
[2]) as below: 

k = (30) 



2,-- QL 



' Qo 2ro ' 



i?2 r^^^ 2ro' 



(31) 
(32) 



7 



Class lb : a = 1 , 6 < 



Q{r) = ±\/^kr^ sec[-v/^V'(r)] 

8TTp{r) = ^ + bk'^r'^ sec^[\^ip{r)] 
The solution set Class lb, are confirmed for 

b<0 ^ M> \Qo\. 

The corresponding coefficients are: 

»-'-(^'' 



.M 



2M QLi,M go^-l.2M 



I = [(tt)' - 1]"^ sec-n(— - ^){{j^r - 1}-^] + (1 - — + ^)n7r - ^r'i^ 



' Qo ro Qo 



Qo ro 



Qo 



3.2 The case of Class II 

In this case Guilfoyl assumed that the metric potential oj^ has the following form: 

—au2 = c + ip(r). 

The interior solutions for this case are: 



a;2 = ^[ln(-a2V^(r))]2 
[ 87rp(r) = |.-^[-ln(-a2^(r))]-i 



And also the constant coefficients: 



Qo ro 



A; = l(M-^)exp[-2( 
rg ro 



ro 



M Qo,-2n 2M 



Qo ro 



7 _ l^-M _ <3()\/1 _ 2M I Qgar/M _ (9o\-l _ /M _ Qo\M _ 2M , Qo\-ll 
*-2lQo ro A-L ro '^'if' ^^Qo 2ro ) ^Qo ro ro ^ ^ ) J 

•-P[-2(f -^)-^(l-^ + f)]- 
We use these interior solutions in Section 5, to form different types of relativistic stars. 



4 Lagrangian formalism for a charged particle moving in RN spacetime 

Now let us consider a Lagrangian like: 

C = ^muiu' -V{r), (42) 

where = ^ is the velocity 4-vector for the test-particle and V{r) is the potential in[TOj One can 
write down this equation for the RN spacetime as: 

£ = ^mguiu'f - V{r), (43) 

in which the diagonal components of the metric ([2]) are used. Therefore the Lagrangian is a function 
like: 

C^Cit,r,9,(t>,i,r,9,^). (44) 

Here the dot stands for differentiation with respect to affine parameter A in geodesic motion. Using 
the RN metric ([2]) yields: 



+ 1 m r2 (sin (e)f - ^ - J (l - 2 ^ + ^) (m^ + ^ 



(45) 



Since we take the geometrical units, we have dX = dr , where r is the proper time. The action in this 
space time, therefore is defined by [22] : 

S = J C{t,r,e,(l),i,r,e,<j))dT, (46) 

turning to the following one for equatorial rotations: 

S= [ C{t,r,<p,i,r,<}>)dT. (47) 



Varying this action, we can obtain the Euler-Lagrange equation of motion in RN spacetime: 

For i = 0,1,3, this leads to three equations which together are the equations of equatorial motion 
(0 = ^) for the test-particle, moving in RN geometry. Using the definition of the potential in pU|) . and 
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the Lagrangian ()45p in ()48p we have: 



(r(r))^ 



+ -1 + 2^ 



(r{r))^ 
(r(r)Y 



dr 



t(r) 



d2 



t t) 



0, 



-m(|:r(r)) 2 



X 1 



+ 



90 



r{r) {r{T)y 



■ m 



mo _|_ 2 qo^ \ f d 



-1 



(r(r)) 

ir)y 



(r{r)f 
mo _ 2 



(r(r)) 



90 



(r(r))^ 



X 1 



90 



r{T) (r(r))2 



-mr(r)(A0(r)) 




mo 



90 



(r(r))^ 



(r(r)) 



(r(r))- 
3 



0, 



(49) 



2 r (r) (^<^ (r)) |:r (r) + (r (r))^ |i<^ (r) = 0. 

Using various interior solutions, in the next section we solve these equations of motion numerically, for 
a particle moving on a relativistic star (Figure [T} , and obtain the shape of possible orbits due to different 
available potentials. 



5 Moving charge on a relativistic star: various interior solutions 

Now let us come back to our original source, i.e. the star with radius r = Rq (Figure [T]). We will illustrate 
the orbits of a massive charged particle, moving around this star. The interior distributions in the star, 
obey the three cases of interior solutions discussed in Section 3. Since we have discussed the interior 
solutions for a star with radius vq , for the total charge amount of the star in Figure [1] we have: 

qo = Q{Ro) = Q{ro - Sr), (50) 

in which Q{r) is one of the interior solutions introduced in Section 3. Also for this star's total mass 
we have [3]: 

mo = ^ ^''V2dr(8Mr) + ^) + 
Using Eq. (j24p this can be summarized as: 



where was determined in ()32p . In the rest of this section we will discuss different sources. 
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5.1 For Class la 



In this case the total charge qq can be derived from the general formula (|5Up using Class la solution, 
namely Eqs. ([281132]) : 



(52) 



The total mass mo has the same relation as ()5ip . We should notice that, in our investigations we 
suppose that the final amounts of charge and mass, Qo and M , have the same ratio ( ^ ) in all considered 
cases, to facilitate the comparison of interior solutions. The potential in (jlOp . for an intermediate angular 
momentum, is shown in Figure [3l Falling in this potential, the test-particle can access three types of orbits, 
due to the total energy E (notated in Figure E]). 




Figure 3: The potential for a test-particle moving on a relativistic star, with the interior solution Class la. The 
illustration is plotted for Qo = 0.85 , q = 0.18 , m = 0.2 , L — 0.225 , tq = 2 , 6r = 0.1 ro . The unit of length along 
the coordinate axis is M . 

1) Periodic bound orbits 



This type of orbit, is often called planetary orbit, specially when the shape of trajectories, is elliptical. 
Figure [H shows some possible bound orbits for different values of initial energy E . 

2) Hyperbolic motion 



For particles coming from infinity, hyperbolic motions or escape orbits are available. That is when 
E — V{r) = possesses only one zero. Figure [5] shows some possible types of escape orbits, for which the 
particle comes to the vicinity of the star, having a definite minimum distance, and then repels again to the 
infinity. 

2) Terminating orbits 
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Figure 4: The possible bound orbits (or planetary orbits) for a test-particle moving on a star with interior 
solution Class la, using different initial points of approach and different initial energies: a) r = , E = ei\ 
b) r = 2.7 , E = ti\c) r = 7 ,E = €2; d) r = 3, E = e^; e) r = 4 , E = €4; i) r = 9, E = 0.36 ; g) 
r = 60, E = 10. 
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d) 

Figure 5: The possible escape orbits (or hyperbolic motions) for a test-particle approaching to a star with 
interior solution Class la, using different initial points of approach and different initial energies: a) r = 2.8 , 
E = e3;h) r = 8 , E = 64; c) r = 10 ,E = e2;d) r = 20 , E = 0.9. 
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Having high energies, when the particle moves in the vicinity of the star, terminating orbits are also 
possible. This type of orbit is unstable and terminates when the particle falls on the star's surface (being 
captured). Figure O shows this type of motion for different energies. 




Figure 6: The terminating orbits (or capture) for a test-particle approaching to a star with interior solution 
Class la, using different initial points of approach and different initial energies: a) r = 2.197, E = €4; b) 
r = r„, E = 0.9. 

5.2 For Class lb 

The total charge of the star, using Eqs. (13311360 is: 

X sec f,/f^-l{sec-i([# - ^]^^)^^ 

(53) 

+ /l_2M + fio!(M _ fio)-i(2 A/ _ Qo)-i 

-Qo^l-(ro-<^r)2Qo(2f -^)ro-3(M - <^r\2 f - f )-!}) ■ 

Utilizing the same data which were used to plot the potential in Figure O can lead us to the potential 
for Class lb, shown in Figure [71 For a particle moving on a star with Class lb as the interior solutions, 
also the so-called three types of orbits are available. The periodic bound orbits (are shown in Figure [8]), 
the scape orbit (is shown Figure [9]) and the capture (is shown in Figure [TO]) . 

5.3 For Class II 

The total charge in this case, using Eqs. (138Mip is: 

go = {-rl + 2Mro - Qg)(2Mro - Ql)Ql^{ro - 5rf 

X (2 - XrlQl + XQlMro - XQl + 2XM^r^ + 2YQlr'^ - iYQlMro + 2YQ^^ 

V [ - In r , ''o^ _/_9 Yr'^D'^ 

^[ ^O^'^y {^Mro-Ql){2Mro-Ql){-rl+2Mro-Ql) ^ "^^Wo 

2Qg(-rg+2Mro-Qg) 

+XQlMro - XQ^ + 2XM^rl + 2YQlrl - WQlMrQ + 2YQl}e '^'o-Qg)^ ) 



(54) 
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Figure 7: The potential for a test-particle moving on a relativistic star, with the interior solution Class lb. The 
illustration is plotted for Qq = 0.85 , g = 0.18 , m = 0.2 , L = 0.225 , rg = 2 , dr = OAvq . The unit of length along 
the coordinate axis is M . 



in which 



^ ^ / 2Mro-Qg-rg ^ _ . / 2r3M - - rgQg _ 4^2j^jy^ + 2ro5rQ2 + 25r2Mro - 5r^Ql 



^0 ' V ^0 



The corresponding potential for the same data used before, has been ihustrated in Figure [TTl According 
to this potential, the equation E — V{r) = , possesses only one zero. Therefore, it is expected that the 
particle exhibits only unstable orbits (escape and terminating orbits). However, for medium and high 
energies, when the particle comes from infinity, the periodic bound orbits are also available. All types of 
possible orbits have been indicated in Figure [l2l 



6 Motion of a charged particle on a Majumdar-Papapetrou star 

Our last discussion, belongs to the stars with a special case of interior solutions Class la, corresponding 
to the Weyl-Guilfoyle spherically symmetric metric potential in Eq. (j26p . From Equation (j27p one can 
confirm that for such metric, the conditions a = 1 and 6 = form a star with equal total mass and charge, 
i.e. 

M = \Qo\. 

This condition for a pressure-less charged perfect fluid, is the Majumdar condition [25t I26| . In this 
case, the relation between the metric potential and the electric potential ip{r) is a perfect square: 

uiipf = {c + iff, 

where c is a constant. In this section we generally investigate the potential for a massive charged 
particle moving on a Majumdar-Papapetrou star. In order to do so, we consider the shrunk star (the inner 
circle in Figure [1]), to apply the interior solution. We use also, the Class la solution with a = 1, 6 = from 
Ref. The solution is: 
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a) 



d) 



b) 




e) 



c) 





g) 

Figure 8: The possible bound orbits (or planetary orbits) for a test-particle moving on a star with interior 
solution Class lb, using different initial points of approach and different initial energies: a) r = , E = e^; 
b) r = 2.75, E = e5; c) r = 7.4 ,^ = £5; d) r = 7.5, ^ = €0; e) r = 2.9, ^ = £7; f) r = r„ , E = ej; 
e) r = 11, S = 0.6. 
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'a) ^ b) 

Figure 9: The possible escape orbits (or hyperbolic motions) for a test-particle approaching to a star with 
interior solution Class lb, using different initial points of approach and different initial energies: a) r = 10, 
E = e5;h) r = 10, E = ej . 




Figure 10: The terminating orbits (or capture) for a test-particle approaching to a star with interior 
solution Class lb, using different initial points of approach and different initial energies: a) r = 2.725 , 
E = e5;h) r = 2.729 , E = ej; c) r = 4.5 , E = 0.3. 
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0.18 



0.16 



0.14 



V(r) 



0.12 



0.10 




Figure 11: The potential for a test-particle moving on a relativistic star, with the interior solution Class II. The 
illustration is plotted for Qq — 0.85 , q ~ 0.18 , to = 0.2 , L — 0.225 , tq = 2 , Sr — 0.1 vq . The unit of length along 
the coordinate axis is M . 




a) 




b) 




c) 




d) 

Figure 12: The possible orbits for a test-particle moving on a star with interior solution Class II, using 
different initial points of approach and different initial energies: a, b) the periodic bound orbits for r = 15 , 
E = €s and r = 10, E = 0.3; c) the escape orbit for r = 20, E = e% and d) the terminating orbit for 
r = 8.3, ^ = £8. 
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Class la : a = 1 , 6 = 



(r)] 



(55) 



8'7Tp{r) 



in which the function 'ip{r) was defined in (j29p . Since M = \Qq\ , for this class of solutions, the constant 
coefficients are as follows: 



in which 



1 



^0 
/2 = 



2|QoL^ IQok 



(56) 
(57) 

(58) 



As before we use the default solution of Q{r) . Here the default solution is the negative one. There 
would be two cases: 



I) The case I = +kR^ (Type 1): 



In this case the total charge of the star is: 

qo = Q{ro -5r) = - 

Using Eqs. (|55ti57p we obtain: 

„ _ (-2ro+|Qo|)(ro-5r)3|Qo| 

90 — 13 

!i) 



k{ro — 5ry 



{l-kR'^Jl 



. 



rf,-2\Qo\rf,+r^Ql+4\Qo\r^5r-2roSrQl-2\Qo\Sr2rQ+5r'^Ql | 
3 J- 



-1 



(59) 



(60) 



The total mass also, can be calculated using Eq. ([?T|) . This leads us to conclude: 

mo = \qo\, 

which is predictable. The corresponding potential is illustrated in Figure [T3j The minus part of the 
potential is of no importance, because we previously decided to consider only the positive energies. 

For E = in Figure [T3l the escape orbits and periodic bound orbits are available, which are illustrated 
in Figures [14] and [I5l 

Also for higher energies, the periodic bound orbits are possible (Figure [T6j) . 



II) The case I = -kR^ (Type 2): 
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Figure 13: The potential for a test-particle moving on a Type 1 Majunidar-Papapetrou star. The interior solution 
is Class la considering the first case I — +kR^ . The illustration is plotted for Qq = —1 , q = 0.18, m — 0.2, 
i = 0.124 , rQ — 2 , 6r = 0.1 rQ . The unit of length along the coordinate axis is M . 





Figure 15: The escape orbits for a test-particle moving on a Type 1 Majumdar-Papapetrou star for: r = 40 , 
E = 0. 
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Figure 16: The periodic bound orbits for higher energies for a test-particle moving on a Type 1 Majumdar- 
Papapetrou star taking: r = 18, E = 0.2. 



According to Eq. ([59|) . the total charge is positive, therefore using Eqs. (|55ll57p we have: 



qo = niQ 



(2ro-|Qo|)(ro-3r)3|Qo| 



r-4-2|Qo|r3+rgQg+4|Qo|rggr-2ro5rQg-2|Qo|gr2ro+^r2Q§ 



+ 1 



(61) 



The potential in this case for the same data is plotted in Figure [T71 as the previous case. This potential 
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Figure 17: The potential for a test-particle moving on a Type 2 Majumdar-Papapetrou star. The interior solution is 
Class la considering the first case I = —kR^ . The illustration is plotted for Qo = 1 ^ Q = 0.18 , m = 0.2 , L — 0.124 , 
To = 2 , 6r = 0.1 tq . The unit of length along the coordinate axis is M . 



allows periodic bound orbits and escape orbit for the test-particle (Figure [T8l and [T9|) . 

One can discover that for the default relation for Q{r) in Eq. (I55p . the total charge of the star, has the 
possibility to be negative or positive. This depends somewhat on the mathematical relationship between 
the constant coefficients / and . As we could see, this affects the test-particle trajectories and change 
its possibilities to have various kinds of motions in each case. 
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a) 



Figure 18: The periodic bound orbits for a test-particle moving on a Type 2 Majumdar-Papapetrou star, 
using different initial points of approach: a) r = , E = eg] h) r = 5 , E = 0.1990 . 




Figure 19: The escape orbits for a test-particle moving on a Type 2 Majumdar-Papapetrou star for: 
r = 2.2, E = eio. 
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7 Conclusion and discussion 



Beside the differences between the shapes of orbits on different types of relativistic stars, there are other 
items, which can be used for more analytic comparisons. One of them is the period of stable orbits 
corresponding to the minimum energy Emin , coinciding the minimum of the potential curve at , Prom 
Equation d?]) and ([5]), for example, the period of stable circular orbits in Class la, will be (see appendix 
A): 



Ti. = 2vr L Z^, (62) 

in which 

and Vm is indicated in Figure [3l We can notice that: 

I lb ^ I la ; 

and therefore 

hi > hh > iia- 

Note that, for the current data we have: 

^11 < ^ib < Tia. 

Another important item, useful to be compared between the cases, is the precession of periastron/apastron 
in periodic bound (planetary) orbits. For stable orbits, one can derive this precession as (see appendix B): 



ha. ha 



A0 = — = 27T[-6q'oui, + ^ + 6mo< - 1 - -^]^ , (63) 



which has been derived for known angular momentum, energy and the radius of stable orbits in Class 
la. We can conclude that: 

A(/.|ib < A</.ia. 

In this article, we considered a relativistic charged perfect fluid as a relativistic spherical star, and 
reviewed the Guilfoyl's interior solutions, through comparing the motion of a massive charged particle on 
such a star. We investigated the effects of these interior solutions on the exterior geometry, for various 
types of relativistic stars, by studying the geodesic motions of a massive charged test-particle. We showed 
that, how changes in distributed charge and mass in a certain volume, could affect the dynamical properties 
of the exterior geometry. We also plotted the potentials, and possible orbits, confirming that each case, 
may exhibit peculiar shapes of orbits. We also discussed the special case of a pressure-less Majumdar- 
Papapetrou star, and illustrated the potentials and corresponding types of possible orbits. 



Acknowledgements: I would like to thank Zahra Gh. Moghaddam, for her interest in this work. 
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A Derivation of the period for stable orbits 

In Equation p3|) we substitute: 



Emin = C{rm) + Jg{rm){m^ + ^). (64) 
The Equation ([7]) and ([8]), for the current case can yield: 

dr g{rm) ' 
d^^ll 

To calculate its period we need the rate of changing in the coordinate time t with respect to (p ^ 
following: 

By chain differentiation we have: 

dt dt/dr Emin/girm) 



del) d^/dr qjrl^ ■ 

Therefore for the period of stable orbits in Class la, (from Eq. (I64p ) we have: 



B Derivation of the precession of periastron in planetary orbits 

We will find this precession here in the following way [27] : 
Substituting u = j, in Eq. (fTT]) yields: 

(^) =-qoU +2mQU -'u(l + -^)H Jl~^ + £2 ■ (^^) 

We consider an approximately stable orbit, with deviations on the trajectory. This helps us compare 
the trajectory with the known circular orbits and simplify the calculations to derive the precession. By 
defining: 

Z — U Umi 

in which z is the deviation from circularity. We substitute this in (j67p . ignoring the terms 0{z^) 
because we concern about nearly circular orbits ( z <^ 1 ) . Considering L = £ii^ and E = Emm for the 
stable orbits in Class la, yields: 



+z[^^ - ^QWm - - 2um + Gmoul- '-^^ + '-^] " " (68) 

la 2 2 la 2 ^ 

z^l-Gq^u'i + ^ + 6moM^ - 1 - ^ 



la la 
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Here, we must have a periodic function, to return the initial (p and u after a period. As usual we 
choose z = ci + C2 cos(a(/> + C3) in which ci , C2 and C3 are constants. However, the coefficient a is not 
equal to one, and in Eq. (|68|) it would be the square root of the coefficient of z'^ (see |27]). Therefore: 

" = i-Qq^m + ^ + Qrnoul - 1 - ^] i (69) 

ha ha 

For acj) = 27r , we have a complete orbit, therefore the change in (p from one periastron to the next is: 
A<^ = — = 2^[-Qqlul, + ^ + Qmoul - 1 - . (70) 

" ha ha 
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